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Abstract 

Introduction: Image guidance for intervention is applied for complex and difficult anatomical regions. Nowadays, it is 
typically used in neurosurgery, otolaryngology, orthopedics and dentistry. The application of the image-guided sys- 
tem for soft tissues is challenging due to various deformations caused by respiratory motion, tissue elasticity and peri- 
stalsis. 

Aim: The main task for the presented approach is continuous registration of preoperative computed tomography (CT) 
and patient position in the operating room (OR) without touching the patient and compensation of breathing motion. 
This approach is being developed as a step to image-guided percutaneous liver RF tumor ablation. 
Material and methods: Up to ten integrated radiological markers are placed on the patient's skin before CT scans. 
Then the anatomical model based on CT images is calculated. Point-to-point registration based on the Horn algorithm 
during a few breathing cycles is performed using a videometric tracking system. The transformation which corre- 
sponds to the minimum fiducial registration error (FRE) is found during the registration and it is treated as the initial 
transformation for calculating local deformation field of breathing motion compensation based on the spline 
approach. 

Results: For manual registration of the abdominal phantom, the mean values of target registration error (TRE), fidu- 
cial localization error (FLE) and FRE are all below 4 mm for the rigid transformation and are below 1 mm for the affme 
transformation. For the patient's data they are all below 9 mm and 6 mm, respectively. For the automatic method, dif- 
ferent marker configurations have been evaluated while dividing the respiratory cycle into inhale and exhale. Average 
median values for FRE, TRE rigid estimation and TRE based on spline deformation were 15.56 mm, 0.82 mm and 
7.21 mm respectively. 

Conclusions: In this application two registration methods of abdominal preoperative CT anatomical model and phys- 
ical patient position in OR were presented and compared. The presented approach is being developed as a step to 
image-guided percutaneous liver radiofrequency ablation tumor ablation. Implementation of the automated registra- 
tion method to clinical practice is easier because of shortening of preparation time in OR, no necessity of touching the 
patient, and no dependency on the physician's experience. 

Key words: percutaneous ablation, preoperative anatomical model registration, respiratory motion compensation, 
image-guided intervention. 
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Introduction 

Image guidance for intervention is applied for 
complex and difficult anatomical regions. Nowadays, 
it is typically used in neurosurgery, otolaryngology, 
orthopedics and dentistry [1-4]. The application of 
the image-guided system for soft tissues is challeng- 
ing due to various deformations caused by respirato- 
ry motion, tissue elasticity and peristalsis [5, 6]. 

Most patients with liver tumors are not candi- 
dates for surgical resections, because of unfavorable 
pathological sites of tumor, widespread metastases 
and poor general physical condition. For these pa- 
tients there is still a chance in minimally invasive, 
percutaneous ablations of lesions. The major task of 
these procedures is accurate needle placement. The 
conventional method for guiding needle insertion is 
ultrasound imaging, but small and deeply situated 
lesions are difficult to identify. The main task for the 
presented approach is continuous registration of pre- 
operative computed tomography (CT) and patient po- 
sition in the operating room (OR) without touching 
the patient and compensation of breathing motion. 

Nowadays, this kind of system is at the development 
stage [7]. One application of the Image-Guided Sur- 
gery Toolkit [8] is CT-based liver tumor biopsy [9]. Rigid 
transformation and the Horn algorithm [10] are used 
to register the preoperative model and the patient's 
position. The needle is tracked by an optical tracking 
system. Fiducial markers, entry and target points are 
manually selected in CT images. Target registration 
error based on Fitzpatrick estimation [11] is calculated. 
Banovac et al. [12] use an electromagnetic tracking 
system to traverse needles into the radiopaque tumor 
target as well as point-to-point registration algorithm. 
The average registration error between the CT coordi- 
nate system and electromagnetic coordinate system 
was 1.4 ±0.3 mm) in the phantom and 1.9 ±0.3 mm) in 
the swine. Maier-Hein et al [13, 14] proposed a CT- 
based, optical tracking system and a real-time defor- 
mation model to continuously estimate the position 
of the target in the liver. A spline-based approach was 
used to capture local deformations between the pre- 
operative CT anatomical model and the patient's posi- 
tion in the OR. The proposed approach was tested 
combining internal and external markers in a custom- 
designed respiratory liver motion simulator and in 
swine tissues. Diagonal arrangement of the needle- 
shaped fiducials yielded the best results. The rest of 
the article presents the proposed registration method, 



which is used to compensate respiratory movements. 
The proposed method consists of two steps: rigid and 
non-rigid registration. 

Aim 

The main task for the presented approach is contin- 
uous registration of preoperative CT and patient posi- 
tion in the OR without touching the patient and com- 
pensation of breathing motion. This approach is being 
developed as a step to image-guided percutaneous 
liver radiofrequency ablation (RF) tumor ablation. 

Material and methods 

Registering a preoperative anatomical model and 
the patient's position starts with a rigid registration. In 
our approach, the initial registration is based on radio- 
logical markers. The presented method uses Beekley 
CT Spots skin markers which are attached to the 
patient's abdominal skin before the CT examination 
and then are used for patient registration. The select- 
ed markers were used because they do not produce 
artifacts on CT images and it is possible to pinpoint 
their position. Markers are manually selected: in the CT 
abdominal examination and on the patient's skin. At 
the very beginning (it was a pilot study) the Polaris 
Vicra tracking system was used and markers on the 
patient's skin were manually selected. We used up to 
ten markers to test various configurations of the mark- 
er's positions relative to the target point. West et al [9] 
demonstrated the importance of configuration of 
markers' positions relative to the target point to the 
registration error. The description of anatomical posi- 
tion of the skin markers is presented in Figure 1 and 
Table I. The following configuration of markers against 
the central marker were tested: all, up-down, left-right 
and diagonal. For each configuration, the central mark- 
er was set as a target. To find rigid mapping between 
two Cartesian coordinate systems, data must include 
three or more corresponding non-collinear points. 
Horn [10] proposed a closed form solution based on 
a least-squares formulation. Optimal rotation and 
translation matrices are found using singular value 
decomposition {SVD) of a correlation matrix: 

SVD(C) = Udiag{oi)V, 
n 

where ^ _ , = 1 ^' ^' is the correlation matrix, Si are 
points in the first Cartesian coordinate system, 77 are 
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Figure 1 A-C. Diagram of patient's abdominal marker positions 



Table I. Anatomical position of skin markers for manual registration 



Marker 
number 


Anatomical location 


Grounds 


1 


Xiphoid process 


Point fixed to bony structures in the midline of the body 


2 


Intersection of midclavicular line and right costal margin 


Anatomical point fixed to bony structures 


3 


Intersection of anterior axillary line and right costal margin 


Anatomical point fixed to bony structures 


4 


Intersection of midclavicular line and left costal margin 


Anatomical point fixed to bony structures 


5 


Intersection of anterior axillary line and left costal margin 


Anatomical point fixed to bony structures 


6 


First control point on the abdominal surface 


Registration accuracy validation 


7 


Umbilicus 


Anatomical point not fixed to bony structures 


8 


Second control point on the abdominal surface 


Registration accuracy validation 



Surgery at the Medical University of Warsaw after 
acceptance of the ethics committee of the Medical Uni- 
versity of Warsaw (GE Light Speed 16 Pro CT Scanner; 
Contrast medium: lomeron 400 or Ultravist 370; Intra- 
venous bolus: 10 ml of 0.9% NaCI, 1 ml/kg contrast 
medium and thereafter 30 ml of 0.9% NaCI all with 
6 ml/s infusion speed). In order to assess the results, 
FLE, FRE and TRE were used. Following the Fitzpatrick 
estimation of registration error [10], different markers 
were evaluated as targets. The results of this study are 
presented in the next section. After the pilot study we 
found a few disadvantages of the manual registration: 
the necessity of touching the patient's skin manually, 
prolonging the time of the procedure, and dependence 
of accuracy on the physician's experience. 

Due to the disadvantages we decided to imple- 
ment an automatic continuous registration. The auto- 
mated method in this case means no need for multi- 



points in the second Cartesian coordinate system, a,- 
are non-negative singular values of the correlation 
matrix and U and V are orthonormal matrices, 



Rot = U 



10 0 

0 1 0 V T 

0 OdetiUV 1 )- 



and 

Trans = f- Rot(S), 

where S, Tare average values of the point coordi- 
nates in the first and second coordinate systems, 
respectively. 

The method has been evaluated on imaging data 
and patients with any liver diseases from the Second 
Department of Clinical Radiology for patients on the 
Chair and Department of General, Transplant and Liver 
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pie manual indication of the position of markers on 
the patient. Position of markers is determined not 
once, but several times during the respiratory cycle. 
We used the videometric tracking system of Micron 
Tracker Claron Hx. 40. Claron Hx 40 is stereo camera 
which tracks a unique black-white template with an 
error of less than 0.5 mm after calibration. We put 
the markers on the patient before CT scans but un- 
like in a manual registration all markers were placed 
on one (right) side of the abdomen for visibility of the 
tracking system. We constantly grab the position of 
the marker (a few times a second) and calculate the 
FRE error between the DICOM coordinate system and 
patient (tracker) coordinate system. The resulting 
transform is organized by the error FRE and the trans- 
form which corresponds to the minimum is treated as 
the result of rigid registration. 

Breathing motion compensation 

Breathing motion compensation is necessary to 
select a proper breathing phase for needle insertion 
and estimate the target movement. 

For breathing motion compensation, a local defor- 
mation field based on displacement of tracking skin 
markers is calculated. The set of marker positions (in 
the tracker coordinate system) during tracking time is 
treated as reference points for the local deformation 
field, based on a spline approach [15, 16]. The defor- 
mation equation: 

Fix) = D T x G(x) + Ax + b 

calculates a new position of point x in space, 
based on spline coefficients matrix D and affine trans- 
formation Ax + b. Based on the position of tracking 
skin markers in the specific time and after a time 
interval, unknown spline coefficients of the current 
deformation field are found and then this field can be 
used to estimate the new target position. To evaluate 
the correctness of deformation one marker is treated 
as a target. Due to TRE estimation [11] usually the cen- 
tral marker is treated as a target. The transformation 
obtained for minimum FRE in the rigid registration 
step is used to translate the target from the DICOM 
coordinate system to the tracker coordinate system. 

Selecting proper breathing phase 
for needle insertion 

For tracking breathing, the marker positions are 
continuously measured by the videometric tracking 



system. Time course of marker positions is analyzed. 
Breathing cycles are detected based on local ex- 
tremes. The breath is the process variable over time, 
so there is a need for averaging it [17, 18]. Amplitude 
of respiratory motion is not measured directly, 
but through the respiratory phase. The filtering signal 
derived from the so-called breathing bag with 
a Gaussian filter to unambiguously determine the 
minimum and maximum is applied. A mean of the 
maximum breathing phase is assigned by the formu- 
la (assuming zero breathing phase for minima) [18]: 



tr-tr 



N" 



fr;'" - tr 



where A/ max is the number of the observed breath- 
ing cycles and ff 3 * and tf" are moments of time 
corresponding to the next extremes of marker posi- 
tion, and respiratory phase, for any time tj, is inter- 
polated linearly [18]: 



(jj-tr 
{ tr-t\ 



for tf" < f) < tr; 



f f — f max \ 
Lj L ( 



4-mm 4-m 

w+i ~~ l i 



(l_gmax). f()r < fr. 



Calculating the breathing phase allows breathing 
data to be grouped from many breathing cycles and 
also data to be compared between patients. 

Results 

Manual rigid registration 

The rigid registration results of a preoperative CT 
anatomical model and a physical patient position are 
presented both for the phantom (Table II) and patient 
data (Table III). To check the dependence of the re- 
sults on the operator and at the same time not to 
extend too much time with the patient, the method 
was evaluated by two operators. 

For the manual registration of the abdominal 
phantom, the mean values of target registration error 
TRE, FLE, FRE are all below 4 mm for the rigid trans- 
formation and are below 1 mm for the affine trans- 
formation. For the patient's data they are all below 
9 mm and 6 mm, respectively. The smallest errors are 
obtained for the target position located in the center 
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Table II. Errors of phantom markers registration using the rigid or affine transformation in millimeters (Polaris Vicra) 



Kind of error 








Error value 










Error type 




Rigid registration 






Affine registration 






Operator 1 


Operator 2 


Operator 1 


Operator 2 




Casel 


Case 2 


Casel 


Case 2 


Casel 


Case 2 


Casel 


Case 2 


FRE 


2.08 


1.91 


1.56 


4.33 


0.52 


0.52 


0.87 


1.04 


TRE marker 1 


3.29 


2.77 


2.60 


7.27 


0.52 


0.69 


0.87 


1.04 


TRE marker 5 


1.56 


1.21 


1.04 


2.94 


0.35 


0.52 


0.35 


0.52 


TRE marker 9 


1.73 


2.77 


2.25 


6.58 


0.52 


0.52 


0.87 


1.04 


TRE marker 3 


2.42 


2.77 


1.91 


5.37 


0.52 


0.52 


0.87 


0.87 


TRE marker 8 


4.50 


3.81 


3.12 


7.62 


0.69 


0.69 


1.04 


1.21 


FLE 


2.60 


2.42 


2.08 


5.54 


0.35 


0.52 


0.69 


0.69 



Table III. Errors of patient's markers registration using the rigid or affine transformation for markers config- 
uration: 8-4-5-2-6 in millimeters (using Polaris Vicra) 



Kind of error 






Error value 




Error type 


Rigid registration 


Affine registration 




Operator 1 


Operator 2 


Operator 1 


Operator 2 


FRE 


7.1 


5.0 


4.7 


3.1 


TRE marker 1 


11.4 


7.8 


7.4 


4.8 


TRE marker 5 


7.3 


5.0 


4.7 


3.1 


TRE marker 9 


11.8 


8.1 


7.6 


5.0 


TRE marker 3 


9.0 


6.6 


5.7 


4.0 


TRE marker 8 


5.0 


3.5 


3.3 


2.3 


FLE 


9.2 


6.6 


5.9 


4.0 



of the marker coordinate system (target 5 in Figure 2 
- Table II and target 6 in Figure 2 - Table III). 

Automatic continuous rigid registration 
and respiratory motion compensation 

The results of automatic registration of the preop- 
erative CT anatomical model and a physical patient 
position are presented below. The anatomical loca- 
tion of markers is presented in Table III. Due to the 
field of view for the Micron Tracking system we 
restricted the area of marker placement to one side 
of the patient's body. The estimation of a new target 



position is calculated following the deformation field 
equation. Table IV presents two cases for a breathing 
phase divided into three groups: 

• exhale (0.85 < breathing_phase < 0.15), 

• middle (0.15 < breathing_phase< 0.3 or 0.7 < breath- 
ing_phase < 0.85), 

• inhale (0.3 < breathing_phase < 0.7). 

Discussion 

There are some patients who demonstrate a little 
difference in rigid registration error in the whole 
breathing cycle. The error of target position estima- 
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Figure 2 A, B. Diagram of phantom abdominal 
marker positions 


Table IV. Errors of patient's markers registration 
for three breathing phase groups 


Breathing phase 


Error type 


Median value 






Case 3 


Case 4 


Exhale 


FRE_rigid 


7.04 


4.26 




TRE_spline 


4.36 


1.08 




TRE_rigid 


0.57 


0.80 


Middle 


FRE_rigid 


7.05 


7.05 




TRE_spline 


4.37 


1.57 




TRE_rigid 


1.39 


0.49 


Inhale 


FRE_rigid 


7.1 


7.73 




TRE_spline 


4.44 


1.91 




TRE_rigid 


1.8 


1.09 


Table V. Errors of patient's markers registration 
for median FRE 


Median FRE 
breathing phase 


Error type 


Median value 




Case 3 


Case 4 


Exhale 


FRE_rigid 


6.81 


2.44 




TRE_spline 


4.54 


1.18 




TRE_rigid 


1.41 


0.80 


Inhale 


FRE_rigid 


7.18 


9.76 




TRE_spline 


4.16 


2.63 




TRE_rigid 


1.23 


0.85 



tion for a specific patient and for rigid and spline esti- 
mations is similar for the whole breathing cycle, 
but generally it rises with an increase in respiratory 




phase difference between preoperative data and 
a temporary breathing phase. Median of rigid fiducial 
registration error was used as a criterion to distin- 
guish an inhale and exhale phase [13, 14]. Table V 
presents results based on this criterion (for the same 
patient as in Table IV). 

Table VI presents registration and target position 
estimation errors for different marker configurations. 
Reconfiguring the markers will change registration 
errors. 

Meier et ai [13, 14] report worse performance of 
target position estimation for left-right marker con- 
figuration for studies in swine. The results obtained 
generally do not support this rule. Due to big dif- 
ferences between patients (Table VII) it is neces- 
sary to treat each patient individually. One possible 
solution is the initial use of more markers (e.g. 8 to 
10) during CT examination and then after nodule 
localization in the CT preoperative anatomical mod- 
el it is possible to test different configurations of 
markers before intervention to select the best one. 
The final number of proposed markers is limited by 
the line of sight between markers and the camera, 
which is distorted by an ultrasound head and abla- 
tion needle. 

Implementation of the automated registration 
method to clinical practice is easier because of short- 
ening of preparation time in OR, no necessity of 
touching the patient, and no dependency on the phy- 
sician's experience. In OR liver tumor ablation is usu- 
ally performed in general anesthesia, so breathing 
is more regular than free breathing. As regards 
the image-guided system for percutaneous liver RF 
tumor ablation, the problem is to synchronize a res- 
piratory phase for which the static preoperative CT 
anatomical model was generated. Using a shape con- 
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serving mattress, the observation of marker move- 
ment and real time FRE calculation helps to select the 
best respiratory phase for needle insertion. 

In this application two registration methods of 
abdominal preoperative CT and physical patient 
position in OR were presented and compared. This 
approach is being developed as a step to image- 
guided percutaneous liver RF tumor ablation. The 
proposed equation of deformation is calculated in 
real time. Assuming regular breathing conducted 
by the respirator, recurrence of patient position 
and short time between CT study and procedure, 
the proposed solution could be used to synchro- 
nize the respiratory phase for which the static pre- 
operative CT anatomical model was generated, 
which is one of the most challenging steps in RF 
ablation. 

In the absence of a direct correlation between 
the movements of the abdominal area and move- 
ment of internal organs to get more reliable regis- 
tration results in the next step we want to apply 4D 
magnetic resonance or computed tomography 
imaging and build an internal organ motion model 
to correlate internal and external breathing motion. 

Acknowledgments 

The authors would like to acknowledge finan- 
cial support from the Polish National Centre for Re- 
search and Development (NCBiR), grant no. LIDER/ 
03/47/L-l/NCBiR/2010 "Development of planning 
system and computer aided minimally invasive 
surgery for hepatic and metastatic liver carcinoma 
localization, diagnosis and destruction". The me- 
thod has been evaluated on imaging data and 
patients with any liver diseases from the Second 
Department of Clinical Radiology for patients on 
the Chair and Department of General, Transplant 
and Liver Surgery at the Medical University of War- 



Table VI. Errors of patient's markers registration 
for different marker configuration 



Markers 


Median FRE 


Error type 


Median value 


configu- 
ration 


breathing 
phase 




Case 3 


Case 4 


z U p U U V V 1 1 


LAI Idle 


FRF riairl 
rr\L llgiu 


8.16 


6.67 






TD T c n i n a 
1 KL bpil 1 Ic 


2.28 


1.17 






TPF riairl 

1 kl rigid 


0.12 


0.74 




1 n h i Lo 
1 1 1 1 Idle 


FPF riairl 
rKt llglU 


9.96 


9.82 






i kl spiine 


3.18 


0.94 






TPF riairl 
1 KL 1 IglU 


5.51 


U.65 


") iH ! an n ~2 
Z UldgUlldl 


LAI Idle 


FRF riairl 
rr\L llgiu 


6.58 


1.66 






TPF c n i n a 

i kl spiine 


4.54 


0.80 






TPF riairl 
1 KL llglU 


1.96 


1.30 




1 n h ^ Lo 
1 1 1 1 Idle 


FPF riairl 
rKt llglU 


7.13 


1.86 






TPF c n i ri a 

i kl spune 


2.95 


0.94 






TRF riairl 
I KL llgiu 


1.44 


i n a 
1.94 


1 Ipftricht 
Z Ic 1 LI 1 gl 1 L 


LAI Idle 


FRF riairl 
rKL llgiu 


4.79 


4.07 






TPF c n i ri a 
1 KL bpil 1 Ic 


0.85 


1.36 






TPF riairl 
1 KL llglU 


U.D 


U.ou 




1 nhale 


FRE rigid 


4.93 


4.21 






TRE_spline 


0.71 


1.34 






TREjigid 


0.68 


0.94 


All 


Exhale 


FRE_rigid 


7.09 


8.58 






TRE_spline 


4.47 


10.24 






TREjigid 


0.47 


1.04 




Inhale 


FRE_rigid 


7.18 


9.78 






TRE_spline 


4.58 


9.07 






TREjigid 


1.08 


1.14 



Table VII. Errors of patient's markers registration for median FRE 



Error type 










Median value 














Minimum 






Maximum 






Mean 






Exhale 


Inhale 


Total 


Exhale 


Inhale 


Total 


Exhale 


Inhale 


Total 


FRE_rigid 


2.27 


3.09 


3.07 


32.8 


52 


51.69 


9.37 


16.75 


15.56 


TRE_spline 


1.07 


1.42 


1.39 


6.92 


112.88 


52 


2.97 


12.82 


7.21 


TRE_rigid 


0.34 


0.3 


0.29 


1.78 


1.23 


1.69 


0.85 


2.43 


0.82 
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